A THEORETICAL ANALYSIS OF A 
YAW DAMPER SERVOMECHANISM 



BY 

LAURENCE B. RICHARDSON JR. 



PRINCETON UNIVERSITY 

AERONAUTICAL ENGINEERING LABORATORY 



REPORT NO. 172 






l£D«cr/ 

13. r^. 



U&> 



,d\iate 



Schoo'. 



' PRINCETON UNIVERSITY 
AERONAUTICAL ENGINEERING LABORATORY 



PAGE 

REPORT 



A THEORETICAL AMALTSIS OF THE 
APPLICATION OF A YAW DAMPER SERVOMECHANISM TO 
IMPROVE lateral OSCILLATORY STABILITY 



by 



Laurence B. Richardson Jr. 
Lieutenant 
United States Navy 



Submitted in partial fulfillment of the 
requirements for the degree of Master of 
Science in Engineering from Princeton 
University, 1951. 



(72 



PRINCETON UNIVERSITY 
AERONAUTICAL ENGINEERING LABORATORY 



PAGE 

REPORT 



The purpose of this analysis is to investigate the use of an 
automatic control to damp undesirable lateral oscillations* The type of 
control under consideration is a yaw damper servomechanism with rudder 
control applied in proportion to yawing velocity. The effect of angle 
of attack and gyro axis angle on components of roll and yaw picked up is 
developed, and the optimum angle of tilt of the gyro axis is selected. 

The problems of maintaining equal effectiveness of the control through the 
speed range, eliminating yaw compensation during steady turning or rolling 
maneuvers, critical time lag, and prevention of feedback to the pilot's 
controls are considered. 

The airplane used for the analysis is an operational jet fighter 
of high performance which has been reported to have -a "snaking," or "Dutch 
Roll" problem. The fli^t conditions Investigated are cruising at 35,000 
feet, 50!^ nonnal rated power at 10,000 feet, and power approach at sea 
level. Calculations in these configurations predict that the aircraft has 
am vinsatisfactoiy "Dutch Rolling" tendency, probably in all flight condi- 
tions. Installation of the proposed damper system will enable the airplane 
to satisfactorily meet the lateral-directional oscillation requirements. 

In addition, the stabilizer will cause the aircraift to be spirally stable, 
and will have the effect of gearing the rudder to the ailerons for improved 
entry into rolling or turning maneuvers. 
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SYMBOLS AIID COEFFICIENTS 

a angle of attack measured from the flight path to the airplane 
reference axis, degrees. 

6 angle between the reference axis and the horizontal axis, positive 
when the reference axis is above the horizontal axis, degrees 

angle of the fli^t path to the horizontal, positive in a climb, 
degrees 

y\ angle of attack of the principal longitudinal axis of the airplane, 
positive when the principal axis is above the fli^t path at the 
nose, degrees 

C angle between the reference axis and the principal axis, positive 
when the reference axis is above the principal axis, degrees 

0 angle between the gyro spin axis and the reference axis, positive 
whjsn the reference axis is above the gyro axis, degrees 

k angle between the fli^t path and the gyro axis, positive when the 
* gyro axis is above the flight path, degrees 

V angle of axirauth, radians 

0 angle of bank, radians 

P angle of sideslip, radians (vA) 

V sideslip velocity along the Y axis, feet per second 

V true airspeed, feet per second 

p mass density of air, slugs per cubic feet 
O" density ratio, p/p^ 

q dynamic pressure, pounds per square foot (l/2 pV2) 

W wei^t of airplane in pounds 
m mass of airplane in slugs, W/g 

g acceleration of gravity, feet per second per second 
14 relative density factor, n/pSb 

m^ factor for dimensionalizatlon, 2b ^i/V '• 

Cj^ trim lift coefficient, M cos T /qS 
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M niach number 
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"a 



I 

z 

o 



I 



X 



wing span in feet 

distance from the center of gravity of the airplane to the center 
of pressure of the fUi, feet 

wing area in square feet 

control surface area for yaw danger in square feet 
vertical tail area in square feet 
rudder area in square feet 

distance from the longitudinal body axis to the center of pressure 
of the vertical control surface 

the operator, d/dt 

period of oscillation in seconds 

time for eunplitude of oscillation to damp to one-half of its 
original amplitude in seconds 

cycles for the eimplitude of oscillation to damp to one-half of its 
original amplitude 



control gearing, or gain, 



5^ j 



amplitude ratio of damper, l/K or r/&,i 
airplane amplitude, 

moment of inertia about the principal longitudinal axis in slug 
feet squared 



moment of Inertia about the principal normal axis in slug feet squared 



moment of inertia coefficient about the principal longitudinal axis, 
o 

moment of inertia coefficient about the fli^^t-path axis, 

(!' cos^'\ + I‘ sln^-iq ) 

X z 



moment of inertia coefficient about axis normal to fli^t path, 

(I* cos% + I' ain^-m ) 
z ‘ X ' 

I* moment of inertia coefficient about the principal normal axis. I /qSb 
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product of inertia coefficient with respect to flight path axis and 
axis normal to flight path, -(!' ” I')3in'>\ cos"^ 

Z X 

rolling moment coefficient ( ) 
yawing moment coefficient ( Y ) 



side force coefficient ( 



lateral force 
qS 



) 



r,T angular yawing velocity in radians per second 
p, 0 angular rolling velocity in radians per second 

dCi 
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dp 
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The standard NACA conventions in regard to directions of axes, angles, 

control deflections, forces, moments, and linear and angular notions are used. 

*• 

V/here dotted quantities are written, they indicate time derivatives, ie, * dt^ 
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INTRODUCTION 

The purpose of this report is to conduct a theoretical investigation 
into the use of a yaw damper servomechanism to in^jrove the lateral dynamics 
of an operational high speed jet fij^ter. Calculations, flight tests, and 
pilots reports indicate that this aircraft does not satisfactorily meet the 
requiren^nts of KAVAER SR-119B, "Specifications For Flying Qualities of 
Piloted Airplanes," in regard to the damping of the lateral-directional os- 
cillation, controls free or fixed. 

High performance requirements of pi^sent day aircraft and those of 
the immediate future are leading to designs in which lateral dynamic flight 
stability is becoming a serious problem. "Dutch roll" at low airspeeds, 
particularly on landing approaches, is annoying, fatiguing, and dangerous 
to the pilot. It makes instrument approaches particularly difficult. The 
present day trend toward sweepback and low aspect ratio causes high effective 
dihedral at low speeds and also reduces the damping in roll and increases the 
adverse yaw due to roll. These factors combine to contribute to unsatis- 
factory lateral dynamics. "Snaking," the control free lateral dynamic oscil- 
lation, compromises the effectiveness of the aircraft in carrying out its 
mission as a gunnery or bombing platform. Satisfactory solutions to the 
problems of "Dutch roll," and "snaking" will not only improve the flying qual- 
ities, but will give the designers more freedom in still hi^er performance 
designs. 

Considerable work is being done ty various investigators on the causes 
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of these unsatisfactory lateral oscillations. Aerodynamic balance of the 
rudder is a probable major cause. Reference 2 traced the poor damping qualities 
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of its test airplane (Navy AD) to the rudder hinge moment cheiracteristics . 
Reference 3 points out as a possible cause the effect of fuel sloshing on 
airplane dyneimlcs. Non-linearity of the directional stability and damping- 
in-yaw derivatives is shown to be a possible cause in reference 4, Random 
separation phenomena such as flow at the tail-fuselage Juncture and at the 
tail pipe have been demonstrated to have important effects. Reference 5 
shows that positive inclination of the principal axis can have considerable 
stabilizing influence. Here wing incidence has an Important effect. It has 
been demonstrated theoretically that atmospheric turbulence can initiate and 
maintain lateral oscillations. 

The designer is faced with an imposing problem in selecting the proper 
hinge moment parameters, which will satisfy control force and other require- 
ments, and still not invite snaking. Reference 2 discusses this problem. 
Irreversible power-boosted control systems offer a promising solution. In 
any case, it has become apparent that some means must be found to augment 
aerodynamic damping. The increment needed is considerably in excess of the 
amount that can be provided by an increase of vertical tail area or tail 
length. It can be shown that the tail volume coefficients should be deter- 
mined by the requirements for static sideslip stability, unsymmetrical power, 
and stalling qualities, and that other methods should be sought to add to 
aerodynamic damping. 

Although hinge moment parameters combined with friction and slack in 
the control systems explain conventional rudder snaking in many cases, others 
are on record in which the rudder has been locked or observed not to move. 
Usually the calculated and tested control fixed periods and times to damp 
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agree quite well. The case under consideration in this paper is one of the 
latter group. The manufacturer felt that the lateral oscillations were trace- 
able to non-linearities in the derivatives and to separation effects. It was 
originally thought that fuel sloshing was a contributory cause, but baffling 
was installed and found to have no noticeable effect. In their tests they 
confirmed the effectiveness of a "de-boost" tab, as reported in reference 2, 
Since there are no known aerodynamic fixes' for the poor deunping in the rudder 
fixed condition, an automatic stabilization device provides the most promising 
solution. 

It is the purpose of this analysis to consider the effect of using a 
servo stabilizing control, and to select an optimum setting of gain and spin 
axis angle for this application. Two systems are already commercially devel- 
oped (by Hinneapolis-Honeywell and Boeing for the XB-47, reference 8, and by 
Lear for use in conjunction with the F-5 autooilot, reference 6). The theo- 
retical investigation is here conducted with the use of one of these types of 
sei*voraechanisras in mind. 

Reference 7 considers stabilizing devices sensitive to various quan- 
tities and concludes that the only one which has the capability of satisfac- 
torily meeting the requirements of reference 1 is an autopilot which supplies 
rudder control proportional yawing velocity. This is the type of yaw damper 
referred to above. This reference eind the above mentioned systems neglect 
the problem of the tilt of the gyro axis. The gyro, which is the rate of yaw 
sensing device, will pick up components of roll as well as yaw, and at only 
one angle of attack will be a pure yaw sensing instrument as is customarily 
assumed. This fact can be used advantageously to effectively gear the rudder 
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to the ailerons by making the yawing moment due to roll positive. This elim- 
inates adverse yaw upon entry into turns and rolling maneuvers. A stabilizing 
device can be designed then that will improve the damping in yaw to meet the 
lateral oscillation requirements, making the aircraft spirally stable at the 
same time, and also improve the response to the rolling control. The problem 
of maintaining equal effectiveness of this control with variations in speed 
is also discussed in this analysis. 
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DESCRIPTION OF THE AIRPLANE AND THE DAMPER SYSTEM 

The airplane under consideration is a single seat, midwing, high 
performance jet propelled fighter with a conventional control system. It 
has a two piece rudder actuated by cables from the rudder pedals, not power 
boosted. The physical characteristics, dimensions, and stability derivatives 
were supplied by the raanxifacturer and checked with independent wind tunnel 
data or methods of calculation as specified by the NACA. 

The damper system to be applied in this case would be known as a 
force type as opposed to a power-boosted type. This means that the servo 
drive has to have an output of the same order of magnitude as a human pilot's 
control effort. In the other case, the seirvo drive can operate into an arti- 
ficial feel system, boost, or irreversible control power control system. Of 
course, the time lag of this force type system is not increased by the in- 
clusion of the time lag inherent in hydraulic nidder boost systems. Consid- 
erable range in power size is available to meet varying requirements. Time 
lags for a force type system can be reduced to less than one tenth of a second. 

The proposed danger system is comprised of three principal parts. 

(l) A rate gyro, a single degree of freedom gyro mounted with its spin axis 

in the XZ plane of the aircraft at some angle to the body reference axis. 

It is constrained by a spring control whose restraint is caused to vary to 

compensate for speed changes. The gyro produces a precession displacement 

and a corresponding electrical signal proportional to the angular rotation 

of the aircraft in space. (2) A controller-amplifier section, which receives 

the rate signals from the rate gyro, amplifies this signal and imparts phase 

/ 

lead. It contains a derivative network which will not pass steady signals 
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from the rate gyro. This neWork prevents the damper from opposing the 
pilot in a steady state maneuver such as a turn or roll. Any transient dis- 
turbances in the steady state maneuver will be damped out, however, and the 
beneficial damper contribution upon entry into these maneuvers passes. The 
gain and the phase lead are adjustable to obtain the optimum settings for 
particular applications. The compensating spring restraint feattore could 
be included in this unit instead of at the gyro. A typical unit would also 
house the rate gyro in its chassis. (3) The servo drive supplies the power 
to move the rudder upon signal from the amplifier-controller. If used in an 
irreversible application (which it is not in this case) the feedback from 
the load, or position follow-up would not be required since the artificial 
feel spring systems will recenter the load when the signal vanishes. 

The total weight of this system could be kept to within fifteen 
pounds, and the space required to less than one third of a cubic foot. The 
gyro should be located near the center of gravity so as to escape the effects 
of body bending. In this case body bending should be a negligible. This 
problem is treated in reference 8. 

It is appropriate that the distinction between an autopilot and a 
stabilizing device should be mentioned. The autopilot has authority over 
the aircraft, presumably by the choice of the pilot. He may overpower it or 
control it through autopilot control settings, but the autopilot has direct 
control of the aircraft. A stabilizing device, such as the proposed yaw 
damper, operates in conjunction with the piloti the pilot, however, has con- 
trol and is superimposing his control motions over the darap)er. The damper 
can be thought of in the same sense as the vertical tail in that it is a petrt 
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of the aircraft which is in continual operation contributing stability. It 
is the opinion of the author that it should be an absolute requirement that 
a stabilizer shall have no feedback to the control system of the pilot i the 
pilot should not be conscious of the damper' s operation except that he feels 
that the flying qualities are improved, and that he is able to fly a more 
hi^ly damped and smoother path. This isolation can be easily accomplished 
in an irreversible control system. In the application under consideration 
it poses a difficult problem that can probably be best solved by dividing the 
two piece rudder between the pilot and the damper. That is, using a control 
surface for the damper alone. Ho solution in which the pilot can feel the 
damper in operation throu^ the rudder pedals should be considered satisfactory. 



I 
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LATERAL DYIIAMIC CHARACTERISTICS 

The first step in the problen is to detennine what the calculated 
lateral-directional oscillation characteristics of the aircraft are. 

Three oonfiguratione of the aircraft have been considered. 

I. The maximum cruise configuration, at 35,000 feet, clean, at a '4ach 
number of .68. 

II. The 505& Normal Rated Power configuration at 10,000 feet, clean, at 
a Mach number of .32. 

III. The Povjer Approach condition with flaps and gear down and with power 
for level flight at a reasonable approach speed. 

These three conditions give a sufficient range of altitude and speed 
to reasonably judge the lateral dynamics of the aircraft with and without a 
damper system. Configuration II was selected because pilots reports have 
indicated that at that particular altitude, at ^^ach numbers below .6, the 
Dutch Roll or snaking has been particularly noticeable. Conditions I and 
III represent the extreme opposite ends of the speed and altitude range in 
configurations of operational interest. 

The equations of motion as presented in reference 9 have been used, 
with the exception that the effects of the product of inertia have been in- 
cluded only in condition I. In this condition is negative.. In the other 
two conditions is positive and it is conservative to neglect the product 
of inertia. For this particular aircraft the product of inertia was found 
to have little effect. These calculations are reviewed in the section on 
calculations. The data is presented in Table II and plotted in figure 2. 

All the points lie in the unsatisfactory region. Flight tests have 
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shovm good correlation with these points. The calculations are control 
fixed, indicating ein unsatisfactory Dutch Roll characteristic. The boundary 
between the satisfactory and unsatisfactory regions of figure 1 of reference 
1 is currently under investigation. The calculated points show the aircraft 
to be in a region which has been called tolerable, meaning tiring or un- 
pleasant, but not necessarily dangerous. The calculations show that a con- 
dition exists which would certainly be annoying and detrimental to the 
carrying out of the airplane's mission. 

Whether this condition is of sufficient severity to wsorrant a change 
or modification in already operational units is doubtful. That is a question 
of opinion of the pilots operating the airplane. All indications are that 
future designs are going to need automatic stabilization! analyses of this 
type will be called for. Applications to longitudinal motions are also in 
use. Future airplanes may well use automatic stabilization both laterally 
and longitudinally with the pilot imposing his control efforts over the 



automatic controls 
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CONTRIBUriOi; OF THE YiW DAl-IPER 

Refeiring to figure lb the relation between the body reference 
axis, the gyro spin axis, and the wind or stability axis is shown. Sub- 
scripts R (reference), G (gyro), and V (stability) are used in the following 
development. 

For simplification the time lag of the yaw damper system is first 
assumed to be zero. It will be later shown that the critical time lag of 
the aircraft-stabilizer system is considerably higher than the time lag of 
any mechanism under consideration. 

The equations of motion are written with respect to the wind or 
stability axes. The derivatives contributed by the damper must therefore 
be related to the V axis for inclusion in the equations. The rudder, how- 
ever, introduces yawing and rolling moments about the reference or body axis 
R. Also the yawing and rolling moments so introduced are proportional to 
the rate of yaw as measured with respect to the gyro axis G. 

Then the development of the contributions of the damper is as follows: 




(C )q * C S = C 
' n B r 

r 



r 



K H' 



G 



(Cjj)g » ^ ^ 




r 



but 




and 
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Cj2 • (C^)v “ (C5 )g cos a + (Cj^)g sin a 
substituting into the above equations, 



[ 4 ^ y cos f 


+ ^(y sin 1 


] cos 


C y cos 1 


+ ^ sin 1 


] sin 


C 4 ^ ^ cos 1 

k 


+ 0y sin 1 


] cos 


[ ^ y COS 1 


+ sin 1 


] sin 


seen to have 


the form of 


total 



therefore 



dC . dC , 

df 



dC, 



djT 

dCo 



A Co « 



'X ° d^ 



dF ^ ^ 



A C 



n 



« K cos I [C + C d a] 
r r 



-Ksinf a] 



A Co 



'V 



A C. 






Kcos^ [C> + C a] 

r r 

K sin I [C(7 + C a] 

r r 



since 






I 
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n 






C = C 

n ^ 2V 2V 

I* r 

2V 



etc. 



AC* 

n 

r 



2V b 

^ K cos I C 

D •' 






AC = 
n 

P 



K sin I C 
b J n . 



neglecting the product o. ^ji » which is small compared to C^ 

6 

V r 

The damper contributions are seen to be directly proportional to 

true airspeed. The rudder power also is seen to vary with altitude from 

Table II. 

The rudder power coefficient has a negative sign, and the rolling 

coefficient due to rudder deflection has a positive sign. It can be seen, 

therefore that if f has a negative sign the increments AC and AC„. will 

’ n^ .C4, 

be positive, and the increments AC and ACj negative. These are the 

desirable signs for these contributions by the damper. It is therefore 

concluded that the gyro axis should be so set that it is always below the 

stability axis, ( f a negative angle) as shown in figure Ic. 

The Increments AC/, and . have little effect on the lateral 

dynamic calculations, and sum [Cj^ + C a] is considerably less than 

r r 

[C + Cfl a] I for these reasons the conservative assumption is made in 
r r 

further calculations that only AC and AC , are supplied by the yaw 
damper. 



The increments added by the yaw damper are then assumed to be 
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AC ® K COS I C 



A C 



ni 



K sin I C 



^0 



n 



The optimum settings of K and | = f(a,D ) are discussed in a 
later section. 

2V 

Reference 10 treats the effect of C (~ C ) on lateral dynamic 

n D ‘ n;v 
P 0 

stability. t-Iaking positive was shown to have a beneficial effect in 
shifting the points on a plot such as figure 4 to the satisfactory area. 

As K euid -f are increased, the positive increment of lengthens the 

period at substantially the same time to damp to one-half amplitude, improv- 
ing the cycles to damp to one-half amplitude. This increment in +C must 
be limited, however, since too great an amount leads to a long period un- 
stable oscillation of the roots that were Identified as the rolling converg- 
ence and the spiral node. Those roots first become a long period stable 
oscillation and then as C is further increased the oscillation becomes 



^0 , 

unstable. This shows that the product VK sin | must be limited so as to be 
satisfactory in configurations as contrasting as I and III. 

The reasons for seeking +C , contributions from the damper are to 
assist in correcting the Dutch Roll oscillations and to oppose the adverse 
yaw introduced by the ailerons and the normal negative derivative (adverse) 
of yaw due to roll. By overcoming the adverse yaw in this manner, entry 
into rolls and tunilng flight can be considerably improved. This benefit 
can only be enjoyed on the entry into these maneuvers, since the steady yawing 
velocity signal will be intentionally blocked to keep the yaw damper from 
opposing a steady turn Initiated by the pilot. 




I 
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C is seen to have its greatest influence in the b, term of the 

1 

stability quart ic. As this derivative becomes more positive it drives b^ 

to zero. As this occurs the rolling and spiral modes first become a stable 

oscillation, and finally a long period unstable oscillation. The limiting 

value of - f for the fixed K was then calculated by setting b^^ to 

zero, keeping all derivatives except C and C , the same, and introduc- 

ing these as f(| ). The limiting P determined in condition I is shown 

in figure 3. Since a varies over a small known range, and f * a - P , 

it is possible to estimate the range of effective | . Sea figures Ic, 3. 

^ and , the non-dimensional forms of the deunper contribu- 

n n ^ 

r p 

tlons, were seen to vary with speed. Therefore in order to maintain equal 
effectiveness of the stabilizing control throu^out the speed range, K 
must be made to vary inversely with speed. One convenient way to do this 



is to apply the inverse of q 



1/2 



to the spring restraint of the gyro. 



thereby modifying the electrical signal from it. The Vs will cancel, but 

-l/2 

a varying factor of o“ ' remains that will have to be calibrated out as 
well as possible. A close approximation to the desirable variations ip K 
with the three conditions is shown in Table III (K to make *5). Com- 

paring these values with those shown in Table II (calculated by using a 
q” ' variation) it is seen that the contributions tend to be too large at 
altitude. Figures 2, 3, and 4 confirm this. Applying dynamic pressure from 
a pitot source seems to be the most convenient way of compensating, however. 

A non-linear spring and bellows with a pitot tube source such as are used 
in artificial feel systems would probably be the best solution. This penalty 
in gad ge try must be paid to realize the benefits of the positive yaw due to 
roll introduced by the damper. Positive C must be limited, and the goal 
shoixld be similar performance at all altitudes and speeds. 
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SETTING GIRO AXIS ANGI£ AITD GAIN 



The question of Just how damped, a motion the pilot would like to 
have is still under considerable discussion. There is no question that in 
a system of this sort it would be possible to add enough damping with a 
stabilizer to the aerodynamic damping to critically damp all lateral oscil- 
lations. However, the standards for the perfornance and stability of the 
aircraft-stabilizer servomechainism are not as yet established. See refer- 
ences II and 12. 

This analysis was made on the assumption that the gain should be 
the minimum required to bring the cycles to one-half amplitude to one half. 
There is evidence to show that the motion will then appear critically damped 
to the pilot. Reference 8. 

A Routh‘s discriminate investigation using K and | as variables 
is possible, but exceedingly tedious. Determination of the optimum values 
of these variables by a plot of vs with K as a parameter is 

shown in figure 3. This approach is relatively easy €ind equally inforniative. 

Referring to figure 3, curves for K = 1 and K * 2 in condition 
III are shown. The higher value of K does not give sufficiently better 
results to warrant the hi^er gain settings. It is in this condition (sea 
level at approximately 1.20 that the basic value of the gain, or 

control gearing, is established. At hi^er speeds (or q's) the gain is 
reduced as already explained. 

Selecting the basic value of K = 1 in condition III, the effective 



value in condition I becomes .5, in condition II, .817. The curve for con- 
dition I for this effective K is shown up to the limiting value of P , 
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determined as discussed in the previous section. In this particular case 
there is no limit to the gyro angle in condition III. In other woi^is with 
gains below 2, too much +C^ cannot be introduced. The effect of holding 
K constant and inclining the gyro axis to the point that damper contribu- 
tion is entirely is shown in figure 4. The points of | = 0 and 

6 o ^ 

3 ® -90 on each curve contrast the effects of stabilizers sensitive to 

pure yaw or pure roll. 

Selecting a gyro angle of 30° and K = 1, the range of values of 
Cl/2 ^ variation of angle of attack of fifteen degrees is shown. This 

setting of 0 and K is taken for further investigation since it is the 
maximum practical gyro angle at this K that will not introduce an unstable 
oscillation due to too much +C in condition I. Some lower angle might 
prove more desirable from fli^t teats. This is a question of how much 
positive yaw due to roll is best for entry into rolls or turns. From figures 
8, 9, and 10 it can be seen that the response to the aileron control is 
radically modified. 

In general, a quick procedure for selecting gyro angle and gain is 
as follows: 

(1) Using the single degree of freedom equation in yaw, a close 
estimate of the total amount of damping required to damp to any given 

is easily obtainable. From this, the amount of AC . that is to be added 

n ijy 

by the damper can be estimated. Knowing the rudder power, the approximate 
gain to get this increment of damping can be fovmd. 

(2) From the adverse yaw and yaw due to roll coefficients the total 
amount of +AC ^ desired can be estimated. Using the gain calculated above. 
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the required | to get the increment in +C^ estimated above can be solved 
for, using the expression obtained under the damper contribution section. 

For the gyro angle, L? « a - | . 

The compensating device on the gyro will keep these increments essen- 
tially constant with speed and altitude. To check to see if an instability 
has been introduced, select a high speed condition where the final ratio of 
. is high (usually at high altitude and speed) and calculate the 
term of the stability quaurtic. 

Final checks using the three degree of freedom equations and involv- 
ing a separate solution of the quartic for each point as was done in figure 
3 can then be made for a critical configuration. The assumption of zero 
time lag will have to be kept in mind and later checked. At this point it 
probably would be desirable to begin fli^t testing the system for the final 
selection of these paramenters. 

In order to prevent any feedback through the pilots rudder control 
system, the yaw damper for this application will probably have to have its 
own rudder area. The calculations in this report have been made on the 
assumption that the whole rudder was available for the yaw damper. This 
can be done when the movements of the rudder caused by the damper can be 
masked from the pilot through a boosted or artificial feel control system. 



C to C 



'nj = Sb 

r r 



c = c — c 

% V Sb i/n^ 

r r r 



and 
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The terras involving C/j have been neglected in the develonment 

r 

of the contribution of the yaw damper. Since 



A C 



n 



K cos 



4 ^ 



n 



AC * K sin I C 



the same results can be obtained by keeping the product of KC constant: 

^ S 



r 

or as the area Sjj. for the damper, is reduced, the gain K is increased. 
Note that the vertical location of the area for the damper is a negligible 
effect. The maximum practical value of K obtainable with this type of in- 
stallation (probably 2 or 3) will limit the minimum size of the damper rudder 
area. 



In summary, the optimum gyro angle is chosen as 30 degrees, and the 
optimum gain as 1 at sea level and approach speed. This setting has the 
following results: 

(1) Moves the points shown in figure 2 from the unsatisfactory 

region to the satisfactory damping region where than .5. 

(2) Adds a sufficient amount of AC , in each case to make the 

n 

aircraft spirally stable in all conditions. 

(3) Adds a +C increment of an approximately constant value to 
all conditions, which is considerably larger than the total adverse yaw due 
to roll and the rolling control. 

The beneficial effect of the additional deimping is shown in the 
response data plotted in figures 5, 6, and 7. The improved response to the 
rolling control is shown in the response data plotted in figures 8, 9, and 
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10 . 

The effect of varying P at a constant K for condition III is 
shown in figure 4. 

Figure 3 shows the selected points in the three conditions, and a 
^l/2 ^ plot for constant K for conditions I and III including the 

range of angles of attack. 

All calculations Involving the contributions from the yaw damper 
were made without using the small angle assumptions or neglecting the influ- 
ence of roll due to rudder deflection. The simplifying assumptions made in 
the analysis of the damper contributions were found to be justified. 



I 



1 



I 
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ADV/JITACE OF YAW DAi-IPSR STABILIZATION 



It was pointed out in the Introduction that aerodynamic damping 
should be augmented by some other means than increasing the tail length or 
area. Increased tail length is limited by operational factors such as car- 
rier elevator sizes for naval aircraft. The amount of vertical tail area 
also has a practical liirdtation which precludes the increases in damping 
that are desired. 



The yaw damper has the advantage of inoreasing C without effect- 
or 

ing the value of C . Tail size can then be chosen without consideration 
of the requirement of hi^ damping, and C can be set from consideration 



of static stability, stalling qualities, etc. Dismissing the possibility 
of adding to conventional tail length, the vertical tail area would have to 
be increased by a factor of about four to achieve the same results in damp- 
ing that the addition of this stabilizer will. 

Reference 8 points out that with the XB-47 it was physically impos- 
sible to change the effective dihedral negatively enou^ to correct the 
lateral dynamics because of ground clearance, and that the tall area would 
have to be at least doubled to add the necessary damping. It was explained 
that in general, while increasing tail size increases damping, the gain in 
smooth flying qualities is questionable. While the tail size increase adds 



V 



to the damping, the excitation of the aircraft due to side gusts Increases 
in the same proportion. 

This XB-47 report uses the servomechanisms approach to stability 

analysis as outlined in references 11 and 12. In one polar stability plot 

White shows that the effect of increasing C and C together (as in the 

""r 
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case of Increasing tail area) results in little or no iraproveTnent in damping 
because the resonant frequency of the aircraft has been increased. 

Consider the single degree of freedom equation in yaw. 



+ C H' »-C S 
z n^ ng r 

r 

This equation may be transformed into the standard servomechanisms 
form using the following definitions: 



damping ratio 






-C . 

nip 

21 CO 
Z 0 



n. 



undamped natural frequency co 



rudder effectiveness 



K * 

h c 






n. 



The equation of motion then becomes 



M^+2'5co4' +(U^'+' ■= 



Kf UJ 

^ o 






and the 



damped natural frequency 



CJO^ a 00 






resonant frequency 



O) 



•“o/l - 2? 



For the usual values of 5 , co , co , and 00 may bo considered 

or n 

equal. White, reference 8, considering the single degree of freedom case, 
sets up his criterion of damping as the ratio of dynamic yaw response occur- 



ring at co^ , to the effective static yaw response considering co^** co^. 



9 



i 
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The measure of the damping then becomes inversely proportional to the square 



of the undamped natural frequency, hence inversely proportional to 

2 2 

since 0) ^ 0) ^ C /l • 

r o n . ' 2 

P 




9 



Also from the single degree of freedom equation. 



r .695 ^ 

^1/2 2tr ' 5 

White states that within the limits of engineering accuracy, the 
single degree of freedom equations may be used for calculations involving 
the design of the yaw damper and to obtain * 

Results from single degree of freedom calculations are included in 
Table III. 



4 




<4 £ 
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TRAfJSIENT RESP01J3S INVESTIGATIONS 

In order to graphically present the contrast between the transient 
response of the aircraft with and without the damper system operating, the 
following calculations were made. For condition II the responses * 

, Eind 0/pQ were found for a disturbance in sideslip effecting each 
of the three equations of motion. The function introduced is a unit step 
function in sideslip applied as a step function at t “ 0, and neutralized 
as a step function at t * 1. Condition II was selected because in this 
configuration pilots have frequently reported unpleasant lateral oscilla- 
tions. The forcing function used is an attempt to simulate a distxu*bance 
that might be initiated by turbulence or by a quick displacement in sideslip 
however introduced. 

The airplane in the original condition is seen to have persistent 
oscillations in all three plots (figures 5, 6, and 7). The period is long 
enough so that the average pilot could damp them out with the proper appli- 
cation of controls; however, requiring the pilot to continually oppose these 
motions is tiring and detrimental to carrying out his mission. There is a 
considerable improvement in yawing velocity and sideslip, but the rolling 
characteristic with the damper in operation should be investigated in flight 
test. Too high a ratio of roll to sideslip is particularly unpleasant to 
the pilot. From data available in other applications, oscillations that damp 
to half amplitude in less than half a cycle appear dead beat to the pilot. 

From this data and the rudder hinge moment parameters, the rudder 
pedal force that must be supplied by the servo power section of the damper 
can be estimated. An estimate of the amount of rudder throw required of the 



V 
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damper rudder surface can also be made. The final setting of the stops 
limiting the throw of the damper rudder surface must be determined experi- 
mentally since the control required will depend upon the characteristics of 
the servo system as well as the airplane responses and the frequency range 
to be damped. 

The response to an aileron step function is presented in figures 

V 

8, 9, and 10. Condition III is investigated since it is probably the most 
unfavorable. Adverse yaw has been included in the calculations. The large 
positive yaw due to roll introduced by the yaw damper is seen to drastically 
change the initial response upon entry into the rolling maneuver. The air- 
pleine no longer starts to yaw in the wrong direction for the turn, and the 
roll builds up rapidly. The sideslip builds up at a slightly faster rate 
also, but not rapidly enou^ to offset the good effects in yaw and roll. 

The derivative network in the yaw damper circuit will have an RC circuit 
with a time constant so designed as to only pass the first portion of a second 
of the contribution of the damper so that it will have the effect of assist- 
ing in the rapid build-up of roll without the adverse initial yaw, and still 
not oppose the steady-state maneuver. In this manner the damper should con- 
siderably improve the response to the rolling control. 

From figure 5. (the response in yawing velocity to the disturbance 

in sideslip) the maximum yawing velocity per degree disturbance is seen to 

bo approximately two degrees per second per degree sideslip. Then ten degrees 

of rudder ivould be required for a gain of one and an initial disturbance of 

five degrees, since S - jd^ /d 'i' 4* . 

^ r 

Considering a ton degree throw in either direction as the limit for 
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danper applied rudder deflection, it can be soon that the "\/ith damper" 
curves shovm in the aileron response plots will be cut off by the rudder 
reaching its stops. This cut-off point will be determined by the amount of 
aileron applied since figure 9 is presented per degree deflection. V/hen 
the product of the gain and the yawing velocity exceeds the available full 
deflection, the rudder will be initially held against the stop with a ten 
degree favorable deflection. The RC circuit will then cut out the damper 
rudder contribution as the steady signal is sensed by the gyro. Should 
there be a tendency to oscillate, the servomechanism will damp the variations 
about the steady maneuver out. 

The curves showing the improved performance with the damper oper- 
ating (figures 3, 9, and 10) will then be modified by servo rudder stop 
setting. Startling improvement in the entry into rolling maneuvers should 
still be obtainable with limited rudder throw. 

Two factors should be here reiterated. First, these calculations 
\^ere made for the maximum value of ©rro axis inclination, and hence positive 
C Increment, possible at the selected value of gain. In other words these 
aileron responses are for the maximum gearing of the rudder to the ailerons. 
Some smaller angle of gyro axis inclination mi^t prove more desirable from 
flight test or more exhaustive calculations. Second, the device introduced 

I 

to assure nearly constant contributions from the damper with variations in 
speed and altitude will automatically reduce the gain at hi^er velocities. 
This should prevent destructive rudder motions as it hits its stop after the 
gyro senses a large rolling velocity. 
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DETERMINATION OF CRITICAL TIMS UG 

In all the calculations to this point the assumption has been 
made that the time lag of the damper system is zero. The critical tine lag 
of the aircraft-autopilot combination is now investigated to determine how 
much of a limitation this assumption imposes. 

The frequency response method outlined in reference 14 is applied. 
The critical time lag of the system is defined as the time lag that results 
in a neutrally stable or steady-state oscillation. The steady-state response 
® airplane to a sinusoidal forcing function of unit ampli- 

tude was first calculated. The calculations were made for condition III 
since in this condition the ratio of the damping to the critical damping 
was lowest. The gain of the yaw damper is known, having been set at K*1 
for this condition. The gain, K. is the aijplitude of the control-surface 
oscillation produced by the yaw damper mechanism in response to the oscil- 
lation of the eiirplane yawing velocity. Intersecting K = |lA with the 

c * 

airplane gain or amplitude plot the maximum permissible time lag in the 
system can be determined. The slope of the line throu^ the point where the 
phase angles are equal is a measure of this critical time lag. The condi- 
tion for stability is that at each angular frequency where the phase angles 
are equal, the ratio of the amplitudes of the airplane to the autopilot must 
be less than one. 

From figure 11 the critical time lag of the system chosen is seen 
to be 1.2 seconds. Critical time lags for gains of two and three are shown 
also. For a damper gain of three, the limiting time lag is ,7 considerably 
in excess of that inherent in typical force type systems. The critical lag 
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shoxm in reference 8 is a little over .5 secon'is for K = 1. For servo- 
mechanisms with larger time lags introduced by boost systems, this calcu- 
lation would undoubtedly show a closer margin. In this case, the time lag 
of the system is not a limiting consideration. 



I 



I 



PRINCETON UNIVERSITV 
AERONAUTICAL ENGINEERING LABORATORY 



PAGE 2.6 
REPORT /7Z 



CONCLUSIONS 

As a iresult of the theoretical investigation of the application 
of a yaw damping servomechanism to the aircraft under consideration, the 
following conclusions are reached; 

(1) This plane does not satisfactorily meet the requirements of 
reference 1 in regard to lateral-directional oscillations. Control-fixed 
calculations predict that the airplane will exhibit a "Dutch Roll" in all 
the conditions investigated to a degree that could be classified as objec- 
tionable but not dangerous. These calculations substantiate pilots reports 
of "snaking" or "Dutch Roll." 

(2) By installing the yaw damper system as proposed, the aircraft 
will then meet the requirements of reference 1. The lateral oscillations 
will appear nearly critically daiaped to the pilot, and he will be able to 
carry out the mission of the aircraft more effectively with less effort. 

(3) The damper will introduce enou^ damping in yaw to make the 
aircraft always spirally stable. It will also introduce a constant positive 
yawing moment due to roll which will Improve the turning and rolling per- 
formeince. The system will be so designed as to improve entry into steady 
turning or rolling maneuvers, but not oppose them. 

(4) Damping introduced by a servomechanism which applies rudder 
control proportional to yawing velocity is probably the most promising if 
not the only satisfactory way to augment the aerodynamic damping. This sta- 
bilizer introduces damping without increasing the resonant frequency (or 
the effective static yaw stability), which insvures results not obtainable 

by increasing vertical tail area. 
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(5) For this particular installation the optimum settings of gyro 
spin axis angle and damper gain are L3 * 30° and K = 1. These values have 
been checked in a limiting condition and the amount of positive yaw due to 
roll introduced is within the limit that would cause a long period unstable 
oscillation. Some lower value of U that would still produce a +C , might 
prove more desirable from flight tests of the effect on rolling maneuvers. 

In all cases, the spin axis of the gyro should be kept below the stability 
or wind axis. Figure Ic. 

(6) It should be required in the case of a stabilizing system 
(as opposed to an autopilot) that no feedback to the pilots controls be 
allowed. In this case the most practical solution appears to be the use of 
a separate rudder area for the yaw damper mechanism, probably splitting the 
two piece rudder between the pilot and the damper. A differential linkage 
or gearing, with the pilot and stabilizer as independent inputs with pro- 
vision for balancing out the servo applied hinge moments could be developed. 

(7) Calculations show that the damper system as proposed will 
improve the transient response to a disturbance in sideslip, and also improve 
the response to an aileron step function input. 

(3) The critical time lag of the stabilizer-airplane system is 
considerably in excess of the maximum time lag to be expected of any servo- 
mechanism of the type oroposed. The assumption of zero time lag in the cal- 
culations proved to be reasonable. 

White, in reference 3, has an excellent presentation for the invest- 
igation of the stability of the stabilizer-airplane combination which he calls 
the inverse frequency response method. It is felt that this type of approach 



I 
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should be used after the experimental frequency response data of the actual 
sei*vomechanisni to be used is available. Otherwise results would be of 
academic interest only. The method presented here will enable the designer 
to quickly select a system with>the required gain, time lag, and power out- 
put from the servo drive and install it near the optimum position relative 
to the airplane reference axis. Fli{^t testing tl^ system from that point 
would probably yield the quickest dependable results unless an sinalog com- 
puter were available. 






4 



I 

I 



f 



I 

I 



I 







i 



PRINCETON UNIVERSITY 
AERONAUTICAL ENGINEERING LABORATORY 



PAGE 3 1 
REPORT /72 



. CALCULATIONS 



The lateral equations of motion are given in reference 9. In all 

cases Cy , Cy . and ^ were assumed tb be zero. The equations including 
P I* 

the product of inertia were used only in case I. Otherwise they were used 
as follows: 



then. 



where 




Cy ) Q 



C,0 + m Dt 
u o 





« 0 
« 0 
» 0 



A ^ + bj + bg + bi A + Lq 



“ A = 0 



b . * m I I 
4 o X z 



tu * -(m C/j, + C„ I )I - m I C 
^ o Xp X z ox 

brt * C^ ^ A T + (mCrt 4" C-^ I }C 4"iuIC *"1110^ 0 

^ "*C« C +niC^C '•ihCaC *• C-C^ I 4* C« C « 

1 ‘■"p" 

b Cj(C|i C •• C/i • C ) 
o L' n^^ i y 



The fastest and most accurate method for the solution of the stability 
quartic was to use a combination of the Lin approximation method and the 
Zimmerman plot as explained in references 11 and 13. With the Lin method, 
values could be obtained that would lead to an intersection in two or three 
points with the Zimmerman plot. In general, there were two real roots and a 
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+C was used there were 
a + ip, 



Vs' 

For the plots in figures 3 and 4, where each point represents the 
solution of a quartic, the equations for ^^4 3 g 1 0 used where the 

variables were C and C . determined from the particular value of 
(a, P ) for that point. 

Using the above three equations of motion, a general solution for 
forcing function inputs of A. B, and C respectively was made. The transient 
responses took the following forms: 

^ ®'l^ % 

6lO “ A 

where 

a, * Cm^I 
^ ox 




complex pair, although where a large value of 
two complex pairs. 

lifhere the roots were ^ j.’ ^2 



. - i693 

^ 1/2 a 

3 = 27T 



seconds 

seconds 



AC I + Bm C -C(mC- + C_ I ) 
Uq X o nt o £x( I X 

p 0 0 p 



a, * A(C(, C - C C^.)-BC + CCp,C 

1 ’'0 "p ‘‘0 “f( ^0 ^P 



a ® BC.C — CC.C^ 
o L np L 









PRINCETON UNIVERSITY 
AERONAUTICAL ENGINEERING LABORATORY 



PAGE 33 
REPORT IIZ 



where 



where 



^ ' AI.I 
^2 



X z 




= -A(I C . 


+ I Cn ) - Cm I 


X 


z£0 p X 


* A(Cj^ C 

y 


- C C, ) + B(C.I 


= -BC.C . 


+ CC,Ca • 


L niy 


L 5-y 



m C . )+Cm Off 

° “ “a 



0 

5(t7 



a_A^ + a, A + a 



So = Bm I 
2 o z 



a, « AC« I - B(m C + I ) + Cm Ca . 

JL a p Z O Z ip O ^ Lp 

a “ -A(C^ C - C , ) + B(C„ C , + m C ) 
o ' Xp ^ ° “p 

-C(m^Cj^ + ) 

In the solution for the disturbance in p, became p^ , and A, 

B, and C were C_ , C« f and C 

^P ^p °p 

In the response to the aileron step function, became and 

A, B, and C were 0, , and C • 

■^c 

a a 

A, and hence the roots ^ a T ip are known. 

The solution for the typical lateral transient response takes the 



Alt 



>2^ 



g(t) = J?''-" + CgC'"^ + Ei.®^ 3in(pt + F) 



form 
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By use of the Heaviside Expansion Theorem, 
a 

C « ^ 

° ^o 



+ SbgAj + 2b^X^ + b^) 



Then from end conditions 



E * 



•<°o h * °2 > 



sin F 



tan F 



-P(C^ + + Cj) 



<=5 ■ >^2=2-^ *<=1^ V 



where represents the initial condition of , P , or 0 at t * 0. 

Cj = ^ . An initially relaxed system is assumed. 

To find the frequency response for the critical time lag calculations, 
standard servomechanisms methods were used. It is believed that the best 
method of reducing the data is outlined in reference 8, and is as follows: 

For the airplane. 



1 

3 



3 2 

£uD + a„D + a,D + a 
— ^ c 1 o 

r b^D^ + bjD^ + + b^D + 



Using the general equations developed previously with forcing func- 
tions A, B, and C as C„ sin cut, Ca sin cut , and C sin JOt the a' a 






s 






can be quickly calculated. 

Then , the steady-state frequency response of the airplane, can 
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be expressed 



K. 



5r 












where the substitution D ® lo) is made. 



X, 



-a^co + 



4 V 2 . ^ 
CO - bftCO + b 

2 o 



-SjCO + a^w 



-1 

0 . « tan ^ 
1 



-b^co + bj^co 



0 ^ * tan 



Y 

-1 ^2 



Then for each frequency the magnitude of the airplane gain or ampli- 



tude is 



1 ;r 1 « 




X^ cos 


I’'aI 


§r 


Xg cos 



and the phase angle is 

0 * 01 “ ^2 
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T/IBLE I 

Principal Dimensions 

S * 250 aq. ft. 

b « 35.25 ft. 

’ 7.45 ft. 

= 13.45 ft. 

3^ = 38 sq.ft. 

= 8.54 sq. ft. 

c . « 6.75 ft. 

vt 

* 7.42 sq. ft. 

V7 •» 12,000 pounds plus fuel 

Total control throw for the rudder, 25°+ 1° each way. 
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TABLE II 





Condition I 


Condition II 


Condition III 




CR max clean 


1 

30% NRP clean 


Landing 


h 

1 


35, 000 


10,000 


aea level 


V ‘ 


660 


259 


183 


M 


.679 


.323 


.134 


e ‘ 


.736(10)'^ ' 


1.756(lof 


2.378(10)'^ 


W ! 

1 


13,500 1 


13,500 


12,500 


^ 1 


64.6 


27.1 


20.0 




160.5 


59.2 


40 ;o 


“o 


6.8 


7.38 


7.71 


b/2V ; 


.0267 


1 

.0681 


1 

.0964 


Cc i 


.337 


1.18 


1.20 


oc ; 


2.9 


8.7 


1 7.8 


^ 1 


3.63 


3.63 

, i 


4.07 

i 


^ ! 


- .73 


5.07 ; 


; 3.73 


D . 


30 


30 

1 


! 30 


1 


- 27.1 


-21.3 


j 

1 -22.2 




7778 


7778 


1 7822 




24244 


24244 


1 24921 


Ix 


.0055 


.0149 


1 .0220 

1 


I 7 


.0172 


.0465 


' ,0688 

j 




.00015 


- .0028 


- .0030 


K 


.5 


.817 


1.0 

j 


Op 


- .546 


-.400 


^ .400 


'^In. 


.132 


.174 


.231 


0 ^ 


- .119 


- .100 


- .064 



231 

064 




I 
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TABLE II (continued) 



r" 


.0103 


.0043 


I 

.0023 i 








.082 1 

1 




- .0146 


- .0272 j 


- .0386 i 




00035 


.0119 1 

1 


.0223 


Grip 


+ .0380 


- .0350 


- .0520 


Go. 


- .230 


- .175 


- .125 ‘ 

j 


n 

0 / 1 ^ 


.135 


.128 


.135 i 




- .0575 


_ .0444 


- .0410 1 

! 


^n<: 






.015 I 


1 


- .0062 


- •0119 ! 


! - .0121 


Gn- 


.0010 


- .0024 


i - . 0050 


Grp 


- .671 

j 


- .562 i 

i 1 


i -.3^4 i 

1 




1 




i .115 


P 

( 


2.26 


3.56 1 


1 4.13 i 




1 3.15 


7.93 


; 8.28 i 

i f 


I C'/j 


1.26 


1.94 


1 2.05 1 


1 m 

1 

L. 


j 

spirally stable 


31.5 


15.1 ! 




With Yaw Damper Operative 




A Gnv)/ 


! 

! - .0255 


-.0341 


-.0379 


; ACr^ 


1 ' .0137 


.0129 


.0155 

t 


j ^*'4' 


j - .0315 


-.0450 


-.0500 j 


i 


i' .0147 


.0115 


.0105 ! 


’ p 


2.44 


4.28 


4.73 j 


1 


.323 


1 . 35 


2.02 1 


j ='/z 

1 


.13 

) 

J_ . 


.31 


.43 
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TABLE III 

Single Degree of Freedom 





Condition I 


Condition II 


Condition III 




.0172 


.0465 i 


.0688 


( 

Cn y j 


- .0061 


- .0119 1 


1 -.0121 

\ 




.135 1 


.128 


.135 


CJq , 

1 
1 


2.805 1 


1.661 


1.404 


I 


2.80 


1.658 


1.40 




.064 


.077 


.063 


P 


2.24 


3.78 


4.48 


T-/2 ' 


3.88 


5.66 I 


1 7.90 




1.72 


1.43 


1.76 


K for 1 

-5 i 

( 


.26 


.51 


.76 


C-/i If 
K= 1 


.13 


.304 ■ 


.402 
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